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Every second greater than 10^® antineutrinos radiate to space from Earth, shining like a faint antineu¬ 
trino star. Underground antineutrino detectors have revealed the rapidly decaying fission products 
inside nuclear reactors, verified the long-lived radioactivity inside our planet, and informed sensitive 
experiments for probing fundamental physics. Mapping the anisotropic antineutrino flux and energy 
spectrum advance geoscience by defining the amount and distribution of radioactive power within Earth 
while critically evaluating competing compositional models of the planet. We present the Antineutrino 
Global Map 2015 (AGM2015), an experimentally informed model of Earth’s surface antineutrino fiux 
over the 0 to 11 MeV energy spectrum, along with an assessment of systematic errors. The open source 
AGM2015 provides fundamental predictions for experiments, assists in strategic detector placement 
to determine neutrino mass hierarchy, and aids in identifying undeclared nuclear reactors. We use 
cosmochemically and seismologically informed models of the radiogenic lithosphere/mantle combined 
with the estimated antineutrino fiux, as measured by KamLAND and Borexino, to determine the 
Earth’s total antineutrino luminosity at 3.4^2 2 ^ lO^^t'e/s. We find a dominant fiux of geo-neutrinos, 
predict sub-equal crust and mantle contributions, with ^1% of the total fiux from man-made nuclear 
reactors. 



FIG. 1. AGM2015: A worldwide 0^ flux map combining geoneutrinos from natural and decay in the Earth’s crust 

and mantle as well as manmade reactor-i/e emitted by power reactors worldwide. Flux units are i/e/cm^/s at the Earth’s 
surface. Map includes Ue of all energies. Figure created with MATLAB’^®. 
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INTRODUCTION 

The neutrino was proposed by Wolfgang Pauli in 1930 to explain the continuous energy spectrum of nuclear beta 
rays. By Pauli’s hypothesis the missing energy was carried off by a lamentably “undetectable” particle. Enrico Fermi 
succeeded in formulating a theory for calculating neutrino emission in tandem with a beta ray^. Detecting Pauli’s 
particle required exposing many targets to an intense neutrino source. While working on the Manhattan Project in the 
early 1940s Fermi succeeded in producing a self-sustaining nuclear chain reaction, which by his theory was recognized 
to copiously produce antineutrinos. Antineutrino detection projects were staged near nuclear reactors the following 
decade. In 1955, Raymond Davis, Jr. found that reactor antineutrinos did not transmute chlorine to argon by the 
reaction: ^^Cl (tj;, e“) ^’’Ar^. This result permitted the existence of Pauli’s particle only if neutrinos are distinct from 
antineutrinos. Davis later used the chlorine reaction to detect solar neutrinos using 100,000 gallons of dry-cleaning 
fluid deep in the Homestake Gold Mine. Reactor antineutrinos were ultimately detected in 1956 by Clyde Cowan 
and Fred Reines by recording the transmutation of a free proton by the reaction (t^e, e'*') ^n^d^ This detection 
confirmed the existence of the neutrino and marked the advent of experimental neutrino physics. 

Almost 60 years later neutrino research remains an active and fruitful pursuit in the helds of particle physics, 
astrophysics, and cosmology. In addition to nuclear reactors and the Sun, detected neutrino sources include particle 
accelerators®, the atmosphere®d^ core-collapse supernovae®^®®, the Earth^^^s^ most recently the cosmos^"^. We 
now know that neutrinos and antineutrinos have “flavor” associations with each of the charged leptons (e, /i, t) and 
these associations govern their interactions. Neutrino flavors are linear combinations of neutrino mass eigenstates 
(^' 1 , and 1 ^ 3 ). This quantum mechanical phenomenon, known as neutrino oscillation, changes the probability of 
detecting a neutrino in a given flavor state as a function of energy and distance. Neutrino flavor oscillations along with 
their low cross section provide a glimpse into some of the most obscured astrophysical phenomena in the universe and 
most recently the otherwise inaccessible interior of our planet. Antineutrinos emanating from the interior of our planet 
constrain geochemical models of Earth’s current radiogenic interior. Antineutrino observations of the modern Earth’s 
interior coupled with cosmochemical analysis of chronditic meteorites from the early solar system allow scientists to 
model the geochemical evolution of the Earth across geologic time. 

Northern Hemisphere 




Southern Hemisphere 


FIG. 2. AGM2015 reactor-i/e flux in the 3.00-3.01MeV energy bin (in logspace color). Flux units are i/e/cm^/s at the Earth’s 
surface. Note the visible 612 oscillations at ~100km wavelength. Figure created with MATLAB"*®. 

Recently, the blossoming field of neutrino geoscience, first proposed by Eder^®, has become a reality with 130 
observed geoneutrino interactions^^’®® confirming Kobayashi’s view of the Earth being a “neutrino star”®®. These 
measurements have constrained the radiogenic heating of the Earth along with characterizing the distribution of U 
and Th in the crust and mantle. The development of next generation antineutrino detectors equipped with fast 
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timing (~50ps) multichannel plates^^ coupled with Gd/Li doped scintillator will allow for the imaging of antineutrino 
interactions. The imaging and subsequent reconstruction of antineutrino interactions produce directionality metrics. 
Directionality information can be leveraged for novel geological investigations such as the geo-neutrinographic imaging 
of felsic magma chambers beneath volcanos^®. These exciting geophysical capabilities have significant overlap with 
the non-proliferation community where remote monitoring of antineutrinos emanating from nuclear reactors is being 
seriously considered^®. 

Antineutrino Global Map 2015 (AGM2015) shown in Figure 1 merges geophysical models of the Earth into a unified 
energy dependent map of i>e flux, both natural and manmade, at any point on the Earth’s surface. We provide the 
resultant flux maps freely to the general public in a variety of formats at www.ultralytics.com/agm2015. AGM2015 
aims to provide an opensource infrastructure to easily incorporate future neutrino observations that enhance our 
understanding of Earth’s antineutrino flux and its impact on the geosciences. In this study we first describe the particle 
physics parameters used in propagating antineutrino oscillations across the planet’s surface as shown in Figure 2. A 
detailed description of the incorporation of anthropogenic and geophysical neutrino energy spectrum from 0-11 MeV is 
presented which allows for the four-dimensional generation (latitude, longitude, flux, and energy) of the antineutrino 
map as shown in separate energy bins in Figure 3. A vertically stratified model of the Earth’s density, shown in 
Figure 4, based on seismological derived density models are combined with a cosmochemical elemental abundances to 
determine the geological signal of antineutrinos. This signal is then constrained by geo-neutrino measurements from 
KamLAND and Borexino and first order uncertainties associated with AGM map are then presented. 


NEUTRINO OSCILLATIONS 


AGM2015 incorporates the known 3-flavor oscillation behavior of antineutrinos. 
3-flavor Pontecorvo Maki Nakagawa Sakata (PMNS) matrix U: 
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where c^- = cos(dij) and Sij = sin(dy ), and 9ij denotes the neutrino oscillation angle from flavor i to flavor j in 
radians. In this paper we assume 


012 = 0.58?;° 

013 = 0 . 152 ;°:™^ 


( 2 ) 


per a global fit by Fogli et al.^^ in the case of 0i2, and by measurements at the Day Bay experiment®^ in the case of 
013. Phase factors ai and a 2 are nonzero only if neutrinos are Majorana particles (i.e. if neutrinos and antineutrinos 
are their own antiparticles), and have no influence on the oscillation survival probabilities, only on the rate of possible 
neutrino-less double beta decay. We assume ai = a 2 = 0 in this work. We likewise assume phase factor 5 = 0, though 
this assumption may change in the future if evidence is found to support neutrino oscillations violating charge parity 
(CP) symmetry. 

The probability of a neutrino originally of flavor a later being observed as flavor /3 is: 
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FIG. 3. AGM2015 Ve flux (i/e/cm^/s/keV) displayed at 6 select energy bins out of the 1100 total AGM2015 energy bins, which 
uniformly span the OMeV - llMeV he energy range. Each energy bin is lOkeV wide. In conjunction with 720 longitude bins 
and 360 latitude bins, the highest resolution AGM2015 map is a 360 x 720 x 1100 3D matrix comprising ~300 million elements 
total. Figure created with MATLAB"*®. 


where E is the neutrino energy in GeV, L is the distance from its source the neutrino has traveled in km, and the 
delta-mass term = mf — rrij, in eV^. The last approximation assumes no charge parity (CP) violation (d = 0), 

causing the imaginary terms to fall out. The * symbol denotes a complex conjugate, and Uij denotes the element of 
the PMNS matrix U occupying the row and column. Equation 3 can be employed to determine the ‘survival 
probability’ of a he of energy E GeV later being observed as the same flavor a distance L km from its source. In 
particular, the Pe-s-e survival probability of most interest to this paper can be expressed as Equation 4: 
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For simplicity we ignore the Mikheyev-Smirnov-Wolfenstein (MSW) effect^^ on neutrinos as they travel through the 
Earth. We use the neutrino mixing angles and mass constants from Fogli et al. 2012^° and Daya Bay 2014^^ to 
evaluate Equation 4 for all source-observer ranges and energies used in AGM, giving us the survival probability of 
seeing each source from each point in the map at each energy level. This is not a trivial task, requiring > 1 x 10^® 
evaluations of Equation 4 for a full AGM2015 rendering. This is broken down into ~ 1 x 10® point sources, ~ 1 x 10® 
locations on the map at which the flux is evaluated, and ~ 1 x 10® energy bins spanning the 0 < < llMeV energy 

range as shown in Figure 3. Equation 4 can best be visualized in Figure 2, which shows the 0i2 ripples in the 3MeV 
worldwide reactor-i>e flux. 


REACTOR ANTINEUTRINOS 

Reactor experiments have proven the viability of unobtrusive reactor monitoring and continue to contribute 
important information on neutrino properties including the possibility of additional light “sterile” neutrino flavors^®. 

We use the International Atomic Energy Agency’s (IAEA) Power Reactor Information System (PRIS)^^ to identify 
and locate 435 known man-made reactor cores in operation at the time of this writing. PRIS categorizes reactors into 
four states: 

• Operational 

• Under Gonstruction 

• Temporary Shutdown 

• Permanent Shutdown 

AGM2015 includes all “Operational” or “Temporary Shutdown” reactors, including many reactors in Japan affected 
by the Fukushima-Daiichi disaster, which are classified as “Temporary Shutdown” rather than “Permanent Shutdown.” 

PRIS shows 435 Operational and Temporary Shutdown reactor cores distributed among 193 sites, with 870 GWth 
total output after load factor considerations, and 72 reactor cores among 42 sites (total 156 GWth at 100% load factor) 
currently Under Construction. The PRIS database reports thermal capacity directly, which is typically about three 
times electrical capacity (most reactors are about 30% efficient in converting heat into electricity). Historical “load 
factors” of each core are used to convert the total thermal capacity to projected current and future thermal power 
output. Load factors account for the down-time related to maintenance and other outages, allowing for AGM20I5 to 
be a reliable estimate of worldwide antineutrino flux in its year of release. Typical PRIS load factors range from 70% 
to 90%. The 3MeV AGM2015 reactor-Pg flux map due to these 435 Operational and Temporary Shutdown reactor 
cores is shown in Figure 2. 

The reactor-i>e energy spectrum assumes the shape of an exponential falloff in energy, with most reactor-Pg released 
at the lowest energies. To obtain an T'g/s/GWth reactor scaling we assume: 

• The mean energy released per fission Ef is around 205 MeV^®. 

• The reactor thermal power Pth is related to the number of fissions per second A(/^ = 6.24xI0^^s“^Pth[GW]/P/[MeV]^®. 

• 6 Pe created per fission®® 

• 2 z/g created per fission on average above the inverse beta decay (IBD) energy threshold of E^, > 1.8 MeV®®. 

These assumptions yield 1.83 x 10®® z?g/s/GWth (of all energies) emitted from a reactor, with 0.61 x 10®° T'g/s/GWth 
emitted above the 1.8 MeV IBD detection threshold. This is a mean value representative of a typical pressurized water 
reactor at the beginning of its fuel cycle®®. We find the summed worldwide reactor power output to be 870GWth, and 
the resultant Pg luminosity to be 1.6 x 10®®Pe/s and 0.04TW. Jocher et is recommended for a deeper discussion 
of reactor-Pg detection via IBD detectors. 


GEONEUTRINOS 

Observations from geology, geophysics, geochemistry, and meteoritics allow for a range of non-unique solutions for 
the composition of the Earth. The relative proportion of Fe, O, Mg, and Si in chondritic meteorites individually 
varies by ~ 15% each and reflects spatial and temporal differences in where these rocks formed in the early solar 
nebula. Likewise, refractory elements have 25% variation in their relative abundance, which translates into a factor 
of two in absolute concentration difference of these elements. Even greater enrichment factors of these elements 
occur when the volatile inventory (e.g., H 2 O, GO 2 , N 2 ) is mostly lost, as during terrestrial planet assembly. Finally, 
because the Earth’s core is taken to have negligible amounts of Th, U and K®®’^®®, due to their limited solubility in 
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CRUST1.0 layer 1: water luminosity 




CRUST1.0 layer 2: ice luminosity xIO^^ 

4 


xIO' 


xIO 
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FIG. 4. AGM2015 luminosity per km^ per GRUSTl.O layer. Each GRUSTl.O layer is composed of 180x360 1° tiles, each 
with a defined thickness (ranging from 0-78 km) and density (ranging from 0.9-3.4 g/cm®). Note the layers have different 
colorbar scales. Figure created with MATLAB^®. 


core-forming metallic liquids, this becomes another 50% enrichment factor in the radiogenic elements in the silicate 
Earth. Consequently, compositional models predict between 10 and 30 ng/g U (and Th/U = 3.9, the chondritic ratio) 
for the silicate Earth. Given the planetary ratio of Th/U and K/U (1.4 x 10^)^°, and the absolute U content of the 
silicate Earth, its heat production for a 10 ng/g U model roughly corresponds to a surface heat flow of 10 TW and 
likewise 30 ng/g U to ^ 30 TW. Estimates of the Earth’s radiogenic heat production thus vary from low power models 
(10-15 TW of power from K, Th, and U), through medium power models (17-22 TW), and to high power models (> 
25 TW)^^. Accordingly, detecting the Earth’s flux of geoneutrinos can provide crucial data to test competing theories 
of the bulk Earth. 

Two observatories, one in Japan (KamLAND) and one in Italy (Borexino), are making ongoing measurements 
of the surface flux of geoneutrinos at energies above the IBD threshold energy > 1.8 MeV. At Japan the flux 
measurement is (3.4±0.8) x 10® cm“^ while at Italy the flux measurement is (4.3±1.3) x 10® cm“^ Note 
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TABLE I. AGM2015 reactor-i/e parameters and oscillation parameters. @12 oscillation parameters from Fogli et al?^ , and 
013 oscillation parameters from results of the Daya Bay"^ Ve experiment. 
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TABLE 11. AGM2015 distribution and properties of U, Th, and K, which are the main emitters of electron antineutrinos. 
Abundances in the various stratified crustal layers are shown at the top of the table, including Oceanic Grust (OC), Continental 
Crust (CC), and Lithospheric Mantle (LM). Relevant isotopic properties are presented at the bottom of the table. Note the 
abundances are unit-less fractions. Uncertainties shown in this table derived from Huang et Arevalo et Gando et 

al}^, and Bellini et al}^. 


that it is sometimes convenient to express geoneutrino flux as a rate of recorded interactions in a perfect detector 
with a given exposure using the Terrestrial Neutrino Unit (TNU)^", however in this work we focus on simple De flux 
/s) and luminosity (7'e/s). 

AGM models the Earth as a 3D point cloud consisting of roughly 1 million points. National Oceanic and Atmospheric 
Administration (NOAA) Earth TOPOgraphical f (ETOPOl) “ice” data^^ is used to provide worldwide elevations 
with respect to the World Geodetic System 84 (WGS84) ellipsoid. Zero-tide ocean surface corrections to the WGS84 
ellipsoid were obtained from the National Geospatial-Intelligence Agency (NGA) Earth Gravitational Model 2008^^ 
(EGM2008) for modeling the ocean surface elevations around the world. Underneath these surface elevations we model 
8 separate crust layers using GRUST 1.0^®, shown in Figure 4, as well as a 9th adjoining layer per Huang et al.^^ 
which reaches down to the spherical mantle, creating a seamless earth model. Gertain crust tiles which are too large 
(about 200 km across at the equator) to be adequately modeled as point sources are instead modeled as collections of 
smaller tiles using numerical integration, which recursively subdivides large tiles into progressively smaller sub-tiles 
until the contribution of each is less than 0.001 TNU. 

Geoneutrino flux is produced from the decay of naturally occurring radioisotopes in the mantle and crust: ^^®U, 
232Th, 235 u^ 40^^ 1130^^^ i38La, i^^Lu, and is^Re^^. However, we only consider and in our 

flux maps as all other elements’ energy spectrum is considerably below the IBD energy threshold of > 1.8 MeV. 
All abundances for the crust and mantle can be seen in Table H. As shown in Table HI, K is the largest contributor 
to De luminosity but its energy is below the IBD threshold. All elements other than U, Th, and K have a negligible 
contribution to the Earth’s Pg luminosity. 

Successful detection of Pg below 1.8 MeV remains elusive; if successful the incorporation of the remaining radioiso- 
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TABLE III. Contribution of geoneutrino luminosities L in AGM2015 for ^®^Tli, and emitted by the Earth. The 

reactor-i/e luminosity is 1.6 ± 0.1 x s and 40 GW for the world’s 435 reactor cores, which together output 870GWth^^. 

Uncertainties shown in this table derived from Huang et Arevalo et Gando et al}^ , and Bellini et al}'^ . 
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TABLE IV. Contribution of geoneutrino luminosities L in AGM2015 above the IBD threshold > 1.8 MeV for ^®^Th, 
and Ve emitted by the Earth. The reactor-i/e luminosity > 1.8 MeV is 0.6 ± 0.1 x 10^®i2e/s and 26 GW for the world’s 435 
reactor cores, which together output 870 GWth^^- Uncertainties shown in this table derived from Huang et Arevalo et 
Gando et al}^, and Bellini et al}^. 


topes would be beneficial to future versions of AGM. The Earth’s core was assumed to have no significant contribution 
to the flux due to limiting evidence for a georeactor^® and no appreciable amount of ^^^Th, or isotopes^^. 
While certain core models support upper limits of K content at the ^100 ppm leveP®, which would be sufficient for up 
to ^1-2TW of radiogenic heating in the present day, “constraints on K content are very weak”^®, and in the absence 
of stronger evidence we’ve chosen to assume a K-free core. 

Mantle abundances were derived from empirical geo-neutrino measurements at KamLAND^^ and Borexino^^. We 
deconstructed the reported geo-neutrino flux from each observation into separate contributions from U and Th ac¬ 
cording to a Th/U ratio of 3.9. From each of these, we subtracted the predicted crust flux contributions^® at each 
observatory, averaging the asymmetric non-gaussian errors, to arrive at estimates of the mantle contributions. We 
then combined the estimates of the mantle U flux and the mantle Th flux contributions from each observation in a 
weighted average. The resulting best estimates for the mantle U and Th flux contributions were finally converted 
to homogeneously distributed mantle abundances using the spherically symmetric density profile of the Preliminary 
Reference Earth Model (PREM)®® along with a corresponding correction to account for neutrino oscillations. Corre¬ 
sponding values for K were found by applying a K/U ratio of 13, 800 ± 1300®®. The resulting AGM20I5 U, Th and K 
mantle abundances are presented in Table II. The main sources of uncertainty in these estimates are the observational 
errors in the flux measurements and limited knowledge of the subtracted crust fluxes. A detailed description of the 
methods and relevant conversion factors used here are presented in Dye®®. 

AGM2015 neutrino luminosities are for total numbers of neutrinos. Although almost all are originally emitted as 
electron antineutrinos, on average only ~0.55 of the total remain so due to neutrino oscillations. We calculate the 
total Earth I/e luminosity to be 3.4?^?'2 ^ s“®. A detailed breakdown of ®®®U, ®®®Th, and geoneutrino 

luminosity from the lithosphere and mantle can be seen in Table III (for all energies), as well as in Table IV for 
Ev > 1.8 MeV. Figure 5 shows the combined AGM2015 crust±mantle Pe flux. 

A comparison of panels in Figure 4 reveals a more luminous upper crustal signature for the Tibetan plateau, 
Himalayan and Zagros Mountains, and the North American Cordillera. Whereas the lower crust of the central Andes 
Mountains is more luminous compared to the aforementioned bright regions of the upper crust which on average are 
less luminous. One potential explanation for contrasting observations might be repeated upper crustal thrusting in 
the Himalayas and like regions, while magmatic additions to the lower crust of the Andean arc could account for their 
relative brightness. These hndings deserve further investigation to understand better their origin and distribution, 
particularly from seismic and geological perspectives. 
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FIG. 5. AGM2015 geoneutrino flux due to and decay in the Earth’s crust and mantle. Flux units are j/e/cm^/s at 

the Earth’s surface. Map includes of all energies. Figure created with MATLAB^®. 


UNCERTAINTY 


The underlying interior structure and composition of the Earth is, in some regards, still poorly understood. The 
concentration and distribution of radioisotopes, whose decay chains produce geoneutrino flux, dominate the uncer¬ 
tainties. Therefore modeling of the distribution, energy spectra, and total flux of geoneutrinos remains a challenging 
task on its own. A full description of the uncertainty in each element of the AGM flux maps is not available at this 
time, however, we have defined the uncertainties in the specific building blocks of the AGM in Tables I and II as well 
as the systematic uncertainties present in the various geoneutrino luminosity categories in Tables III and IV. The 
uncertainties in Tables I and II in particular can be used to create Monte Garlo instances of the AGM flux maps, 
which could be used to evaluate the variances in each element of the AGM flux map, as well as co-variances between 
map elements. Such a full-scale Monte Garlo covariance matrix is impractical, however, due to the large number of 
map elements, ~300 Million, which would end up producing a 300 Million x 300 Million full-size covariance matrix. 

In most parts of the AGM2015 map uncertainties are strongly correlated over space and energy. This is due to the 
fact that the greatest uncertainty lies in ingredients that affect all map elements nearly equally, such as elemental 
abundances in the Earth’s crust and mantle. Other ingredients which might introduce more independent uncertainty, 
such as the volume or density of specihc crust tiles, are much more likely to introduce minimal correlations, or only 
slight regional correlations. Since reactor-i/e flux is generally better predicted than geo-P^ flux, and since near a core 
the reactor-i?e flrrx will dominate the overall flux, we would expect smaller fractional uncertainties near reactors 
than in other regions of the world. 

We attempt to apply systematic uncertainties to AGM2015 where appropriate, and to derive these uncertainties 
from previous work in the field where possible rather than reinvent the wheel. Eor our uncertainty models we turn to 
Huang et al.^^ and Dye^^. Doing so allows us to apply systematic uncertainties to the various source categories in 
Tables III and IV while avoiding the significant computational burden of a full Monte Carlo analysis, as well as the 
questions that would arise afterward of how to describe the various levels of regional correlation in flux uncertainty 
across space and energy. 
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CONCLUSION 

Electron antineutrino measurements have allowed for the direct assessment of 7-29TW^^ power from U and Th along 
with constraining a geo-reactor < 3.7 TW at the 95% confidence level^^. Such measurements promise the fine-tuning 
of BSE abundances and the distribution of heat-producing elements within the crust and mantle. Several such models 
of the Earth’s antineutrino flux^^do,4i existed before the observation of geoneutrinos; with several recent models being 
presented with the inclusion of geoneutrinos^®’^^d3 Qf aforementioned models incorporate several geophysical 
models based on the crust and mantle from traditional geophysical measurements (seismology, chondritic meteorites, 
etc.) This effort, AGM2015, aims to consolidate all these models into a user-friendly interactive map, freely available 
to the general public and easily accessible to anyone with a simple web browser at www.ultralytics.com/agm2015. 

Euture work includes completing a more detailed uncertainty study using Monte Carlo methods. Such a study 
requires an accurate understanding of the uncertainty in each of the ACM elements listed in Tables I and II, however, 
as well as the correlations that govern their interactions. Future geo-Dg measurements, updated flavor oscillation 
parameters, advances in crust/mantle models, and the ongoing construction and decommissioning of nuclear reactors 
around the world necessitates a dynamic ACM map capable of changing with the times. For this reason we envision 
the release of periodic updates to the original AGM2015, which will be labeled accordingly by the year of their release 
(i.e. “AGM2020”). 


METHODS 

AGM2015 uses CRUSTl.O^® to model the Earth’s crustal density and volume profile via eight stratified layers. 
Elemental abundances for U, Th and K, and isotopic abundances for ^^^Th and for each layer were defined 
by Huang et al.^^ . These values were coupled to well known isotope half-lives and multiplicities to create Vf. luminosities 
emanating from each crust tile. A similar approach was taken with the Earth’s mantle, with elemental abundances 
derived via estimates of geo-7'e flux at KamLAND and Borexino"^^ and density profiles supplied via PREM^®. 

Man-made reactors were modeled via the IAEA PRIS^^ database, with reactor-Pe luminosities found to scale as 
1.83 X 10^° Pe/s/GWth. Reactor-i>e spectra were modeled as exponential falloffs from empirical data^®, while 
®®®Th and spectra were modeled based on the work of Sanshiro Enomoto®^. 

Luminosities from each point-source were converted to fluxes at each map location via the Pe^e survival proba¬ 
bility shown in Equation 4, and a full understanding of the source spectra of each point-source enabled a complete 
reconstruction of the observed energy spectra at each map location. “Smart” integration was applied where necessary 
to more accurately portray crust and mantle tiles as volume-sources rather than point-sources. All modeling and 
visualization was done with MATLAB^®. Google Maps and Google Earth multi-resolution raster pyramids created 
with MapTiler^®. All online content available at www.ultralytics.com/agm20I5. 
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